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ABSTRACT 



ps) ! Turbulence plays a major role in the formation and evolution of molecular 



clouds. The problem is that turbulent velocities are convolved with the den- 
sity of an observed region. To correct for this convolution, we investigate the 
relation between the turbulence spectrum of model clouds, and the statistics of 
their synthetic observations obtained from Principal Component Analysis (PCA). 
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p ' We apply PCA to spectral maps generated from simulated density and velocity 

Q I fields, obtained from hydrodynamic simulations of supersonic turbulence, and 

^ I from fractional Brownian motion fields with varying velocity, density spectra, 

ci '. and density dispersion. We examine the dependence of the slope of the PCA 
pseudo structure function, apcA, on intermittency, on the turbulence velocity 

Kj^ , (/3^) and density (/3„) spectral indexes, and on density dispersion. We find that 

0> I PCA is insensitive to /3„ and to the log-density dispersion as, provided a^ < 



2. For (Ts > 2, a PC A increases with a^ due to the intermittent sampling of the 



"^ I velocity field by the density field. The PCA calibration also depends on inter- 

^ I mittency. We derive a PCA calibration based on fBms with cts < 2 and apply it 

to 367 ^^CO spectral maps of molecular clouds in the Galactic Ring Survey. The 
average slope of the PCA structure function, (apcA) = 0.62 ± 0.2, is consistent 
with the hydrodynamic simulations and leads to a turbulence velocity exponent 
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of (/3^,) = 2.06 ±0.6 for a non- intermittent, low density dispersion flow. Account- 
ing for intermittency and density dispersion, the coincidence between the PCA 
slope of the GRS clouds and the hydrodynamic simulations suggests /3y ~ 1.9, 
consistent with both Burgers and compressible intermittent turbulence. 

Subject headings: ISM: clouds - ISM: kinematics and dynamics - turbulence - 
molecular data 



Introduction 



Turbulence plays a major role throughout the entire lifetime of a molecular cloud, from 
its formation to its fragmentation and collapse into clumps, cores, and stars. Supersonic 
turbule nce is not only respon sible for most of the kinetic energy supporting clouds against 



gravity (jWilliams et a. 



molecular clouds (e. 



2005 



Heitsch et al. 



^ 



.II2OOO ), but is also intrinsically linked to the formation and structure of 



Mac Low &: Klessen 



2005 



2004i Elmegreen fc Scalol 



Vazauez-Semadeni et al.ll2006[ boOTJ 



2004 



Audit fc Hennebelle 



Dobbs et al. 



2006 



Tasker fc Tan 



2009 : iBaneriee et al.ll2009l : Iciover et al.lboiol : Klessen fc Hennebellel I2OI0I : IPederrath et al. 



2OIII ). Turbulent shocks create local, over-dense regions within molecular clouds that may 
collapse into sheets, filaments and pre-stellar cores under the effect of self-gravity, a phe- 
nomenon known as turbulent fragmentation. In this context, the length scale of the en- 
ergy source that drives turbulenc e can explain differe nces in star- formation efficiency and 
type — clustered versus isolated (JKlessen et al.ll2000l ). In addition, the Initial Mass Func- 
tion (IMF), which describes the relative probability of stars with different masses when 
they fo rm, may directly result from the magnitude and spatia l scale of turbulent fluctu- 
ations (jPadoan &: Nordlundl 120021 : iHennebelle fc Chabrierl |2008| . |2009| ). Hence, most ana- 
lytic models of star formation rely on the index of the turbulent energy spectrum (i.e., the 
Fourier spectrum of the velocity fluctuations), and on the Probability Distribution Func 



tion (PDF) of the turbulent density field (e.g.. lElmegreerull997l : IPadoan fc Nordlundl 12002 



Krumholz fc McKeell2005l : IHennebelle fc Chabrien l2008l . |2009| ). The energy spectrum is de- 
fined as i?(k) = 47r/c^ |v(k)| , where v(k) is the Fourier transform of the velocity field, and 
obeys a power law {E{k) oc k'^""), where /3„ = 5/3 for incompressible (Kolmogorov) turbu- 
lence (JFrischlll995l ). and /3j, = 2 for press ureless (i.e., highly supersonic) s hock-dominated 
turbulence, also called Burgers turbulence (lBurgerslll974l : iPassot et al.lll988l ). Intermittency, 
characterized by fiuctuations occurring sporadically (both spatially and temporally) in the 
turbulent fiow, also affects the nature of the turbulent cascade from large to small scales. 
Numerical simulations show that intermittency is mainly caused by the interaction of strong 
shocks, causing rare, strong density enhancements. Manifestations of intermittency are ob- 
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servable in the PDF of the density field ( JKlessenl l2000l : iKritsuk et al.l 120071: ISchmidt et al. 



2009 



Federrath et al.l 12010), and in the index /j,, of the energy spectrum (DubruUd 



1994 



She fc Levequelll994l : iBoldvrevI 120021 : iBoldvrev et al.l 120021 : ISchmidt et al.l 120081 ) . 

While Kolmogorov (incompressible) turbulence does not seem appropriate for molecular 
clouds, in which turbulent motions are highly supersonic (hence compressible), the applica- 
tion of Burgers turbulence is limited to compressive fields for which Vxv = (e.g., expanding 
shells). Three-dimensional numerical simulations of both decaying and driven turbulence, 
however, show that the ratio of co mpressive to solenqidal (for which V • v = 0) v elocity dis- 
persion is 7 = (v^)/(v^) = 0.1-1 JPorter et al.lll998l . Il999l : iFederrath et~aDl201ol ). where Vc 
and Vg are the compressive and solenoidal parts of the velocity field, respectively. An alter- 
native, wliich_jncludes intermittency an d compres s ibility effects, wa s prop os ed on the basis 



(I2OO2I); 



Schmidt et al. 



of the IShe fc Leveaue (1994) model by JDubruUel fll994h : JBoldvrevI fl2002h : JBoldvrev et al. 



(I2OO8I ). In this theoretical context, the turbulent energy cascade in 
the inertial range exhibits properties of Kolmogorov turbulence, while close to the dissipative 
range, intermittent shock structures resulting from compressibility effects start to dominate 
energy transfer mechanisms. Depending on the dimension of the dominant dissipative, in- 
termittent structures, log-Poisson models predict f3y = 1.74-1.83. 

The type of turbulence (Kolmogorov, Burgers, log-Poisson), along with the density 
spectrum and density PDF, infiuence the resulting stel l ar ma ss spectrum, the IMF (e.g.. 



Padoan fc Nordlund 2002; Hennebelle &: Chabriei 



better constrained with extinction data ( Cambresyl 



2008 



1999 



20091). While the density PDF is 



2010a ) , the density spec trum has been investigated both with FIR dust er nission observations 



e^g,. 



1999 



Ossenkopf et al.ll200ll : JBrunt et al. 



Elmegreen et al. 



Block et al.ll2010) and with Hi 21 cm spectral line observations f lStanimirovic et al. 



20011 ). The index of the energy spectrum, A ,, is the subject o 



spectral line studies of mo 



Ossenkopf fc Mac 



low I2OO2 



Brunt 



e cular clouds (JHever fc Schloerb 



2OO3I: JHever fc Bruntll2004J : JHilv-Blant et al.1120081 ). in 



997: Brunt fc Hevei 



many 



2002b 



Heyer fc Schloerbl (119971 ) proposed to use Principal Component Analysis 



eluding this work. 

(PCA) to determine the index of the energy spectrum of molecular clouds using molecu- 
lar line observations. PCA provides pairs of spatial and velocity scales, or "PCA pseudo 
structure functions" , which describe the amount of velocity fiuctuations contained within an 
eddy of a given spatial scale. PCA pseudo structure functions thus provide a description of 
the turbulent energy cascade. Since knowing the turbulence energy spectrum of molecular 
clouds has great implications for our understanding of fragmentation and star formation, 
our ultimate goal is to derive the exponent of the velocity power spectrum from the anal- 
ysis (e.g., via PCA) of spectral line maps. PCA pseudo structure functions, however, are 
based on emission, and hence on a complex convolution between the density, velocity, and 
temperature fields, while a complete description of a turbulent gas fiow requires separate 
statistical descriptions of the velocity and density fields. For instance, turbulent velocity 
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fluctuations can make two physically separated emitting elements at different distance s from 
the observer seem to overlap in velocity space (JBallesteros-Paredes fc Mac Lowll2002l ). This 
would produce an element of double the emissivity, which could be interpreted as a sin- 
gle emitting element of double density or column density. Temperature fluctuations also 
produce emissivity variations that can be interpreted as column density variations. This 
is not only a limitation of PCA, but also of all observational studies aimed at deriving 
the energy spect rum from various mathematical tools, such as Centroid Velocity Incre- 
ment s (see, e.g., Hilv-Blant et al. 20081: Brunt fc Mac Lowl 120041 ). Velocity Channel Anal- 



ysis (Lazarian fc Pogqsvanl l2000t lEsquivel et al.ll2003[ ). and Velocity Coordinate Spectrum 



( iLazarian fc Pogosyanl |2006| ). Unfortunately, while observations provide column densities 
and their PDF, the temperature, velocity, and volumetric density flelds, however, are not 
directly observable. In order to determine the relation between the PCA pseudo structure 
function and the Fourier spectrum of the velocity field, we therefore need to disentangle the 
density, velocity, and temperature contributions to the observed position-position-velocity 
data and to the shape of the PCA pseudo structure function. This can be done by comparing 
the PCA pseudo structure function of simulated spectral maps obtained from simulations 
of molecular clouds to the statistics of the input velocity, density, and temperature fields. 
Here, we use isothermal simulations, and hence neglect the molecular excitation problem and 
the effects of star formation on the gas temperature to concentrate on density and velocity 
fluctuations. Thus, our goal is to establish a calibration relation between the slope of the 
PCA pseudo structure function and the slope of the turbulence velocity spectrum, and to 
inve stigate how this ca libration relation varies with properties of the density fleld. 



Brunt et al.l ( 120031 ) established a calibration relation based on MHD simulations. Their 



simulations included many of the important physical processes in the ISM, such as grav- 
ity, magnetic flelds, star formation feedback, heating and cooling. The focus of their study 
was to link the PCA derived relationship between velocity differences and spatial scale to a 
particular order of structure function. They were able to show that PCA pseudo structure 
functions correspond to a low order (~ 1) structure function even in the regime of strong in- 
termittency. The resolution of their simulation was however too small to allow the existence 
of an inertial range, where the slope of the power-law velocity spectrum can be measured. 
As a result, their simulation had to be modifled a posteriori in Fourier space to create a 
power-law velo city spectrum. In addition to this limitation, the turbulence in the MHD 
simulations by iBrunt et al.l (120031 ) were mostly driven compressively by expanding shells, 
occurring in star-forming regions. 

In this paper, we investigate the statistical relation between the PCA pseudo structure 
function and the intrinsic index of the velocity fleld for two types of forcing — compres- 
sive and solenoidal — based on hydrodynamic simulations with a distinct inertial range 



( JFederrath et al.ll2010[ ). combined with Fractional Brownian Motion (fBm) simulations. We 
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examine the sensitivity of this relation to the statistics of the density field by varying both 
the exponent of the power-law density spectrum and the density PDF of the fBm simulations. 
Section [2] describes PCA and the details of the method. Section [3] describes the simulations. 
In Sections H] and \5[ we present the results of PCA applied to the hydrodynamic and fBm 
simulations, and provide a PCA cahbration in Section [61 We apply P CA to spectral obser- 



vations of molecular clouds taken as part of the Galactic Ring Survey (I Jackson et al.l 120061 ) 
in Section [3 Section |8] consist of a brief conclusion. 



2. Principal Component Analysis (PCA) 
2.1. Method 

PCA, first suggeste d as a tool to derive energy spectra from spectral observations by 



Heyer fc Schloerbl ( 119971 ). detects hne profile differences due to the turbulent nature of the 
flow as a function of spatial scale. The line profile differences are represented by the eigen- 
spectra. The spatial scales over which those line profiles differ are detected in the inte- 
grated intensity images (the principal components) of the eigenspectra. PCA thus provides 
pairs of spatial and velocity scales detected in a PPV cube, similar to a structure function, 
6v = f{Si). This so-called PCA pseudo structure function describes the amount of kinetic 
en ergy contained within a given spatial scale. The details of the method have been described 



m 



Brunt &: Heyerl (l2002al ) and are summarized below. 

The observational data obtained from spectral line mapping of molecular clouds can be 

represented as PPV cubes, which consist of an ensemble of A^,. spectra T(r, v) = T{xi, i/i, Vj) 

= Tij {i = l..Nr, j= l..Ny) of a molecular spectral line at different positions r on the sky. 

Ny is the number of velocity channels. The identifi ed line profile differences are the eigen- 



vector s (or eigenspectra) of the covariance matrix S ( jHeyer fc Schloerblll997l : iBrunt fc Heyer 



2002al ). such that Su*-"'-' = X^"''u^"'\ where n is an index that labels the eigenvectors in order 
of decreasing corresponding eigenvalue (n = 1..^^,). The projection of the eigenvectors u^"^ 
onto the PPV cube ordered by decreasing eigenvalue A*-"-* yields the N^ principal components 
PC(") of the PPV cube, such that PC^") = *Tu(") (of size Nr). 

The spatial and velocity scales characteristic of the principal components and eigen- 
vectors are computed from their Auto- Correlation Function (ACF). The spatial scale 6i at 
which the ACF of the n^^ principal component falls by one e-fold defines the n*'^ character- 
istic spatial scale. In a similar way, the corresponding n*'^ velocity scale Sv is determined 
from the e-fold of the ACF of the n*'^ eigenvector. In the end, PCA provides N^ pairs of 
increasingly smaller spatial (Si) and velocity (Sv) scales, similar to a structure function. This 
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relationship can be approximated by a power law, 5v = voSi"''^'^^ , where vq and apcA are the 
amphtude and the slope of the PCA pseudo structure function. 



2.2. Definition of spatial scales 

The spatial and velocity scales are chosen to be the scales at which the normalized ACF 
falls by one e-fold. In numerous cases, the ACF of the O*'* principal component does not fall 
below one e-fold due to the overall correlation of the bulk of the emission. In such instances, 
the O*'^ principal component was excluded from the PCA structure function. 

Furthermore, the determination of the velocity scales is straightforward because the ACF 
of the eigenvectors is one- dimensional. However, the ACF of the principal component images 
is two-dimensional. Because of the intrinsic shape of molecular clouds, the ACFs of the 1*** 
and 2"*^ principal components are often elliptical, such that the spatial scales depend on the 
direction along which they are calculated. In order to measure spatial scales consistently, 
independent of the geometry of the cloud, the ACFs of the principal components were fitted 
to an ellipse. The ACFs were then rotated such that their long axes are horizontal, and the 
one e-fold spatial scales 6x and 6y were computed along the x and y cardinal directions. The 



spatial scales of the PCA pseudo structure function are then defined as 



to stay c onsistent with pre vious PCA studies (e.g., iHeyer fc Schloerblll997l : iBrunt fc Heyer 



2002a 



bl: lBrunt et al.ll2003h . 



Simulations 



3.1. Hydrodynamic simulations 

We numerically modeled isothermal driven turbulence on a periodic un iform grid with 
1024^ grid cells. We refer the reader to IPederrath etaP J2008I . bood . boioh for the details 
of the simulations. Two kinds of forcing were implemented: solenoidal (or divergence-free) 
forcing for which V ■ f = 0, and compressive (curl-free) forcing for which V x f = 0. These 
two types of forcing mimic actual mechanisms responsible for driving turbulence in the ISM. 
For instance, galactic shear corresponds to solenoidal forcing, while supernova explosions are 
a compressive way to drive turbulence. For each forcing case, 81 snapshots of the velocity 
and density fields, spanning eight large-scale turbulent crossing times were recorded in the 
regime of fully developed, supersonic turbulence. To facilitate the computation of structure 
functions and PCA, all snapshots were resampled on a 256^ grid. The res ampling from 1024^ 
to 256^ did not affect the inertial range scaling (JFederrath et al.l 120091 ). The simulations. 



-7- 



originally with a mean density of unity and an RMS Mach number, Ai = 5.5, were rescaled 
to a velocity standard deviation of 1 km s~^ (T = 10 K, assuming that the gas is composed 
of pure molecular hydrogen) and a mean density of 500 cm~^. 

The statistics of the solenoidally and compressively forced simulatio ns (in particular 



veloci t y and dens ity power spectra and density PDFs) were computed in iFederrath et al. 



(120081 . l2009l . I2OIOI ) and are summarized in Table [H The inertial range of the hydrodynamic 
simulations only extends between k = 5 and A; = 15 due to numerical viscosity. In addition to 
the average statistics summarized in Table (H we computed the inertial range exponent /3„ of 
the energy spectrum for each snapshot. An important point for the upcomin g analysis is that 



the density PDFs of the simulated fields are only approximately lognormal ( IFederrath et al. 



20081 . I2OIOI ) due to significant intermittency. Deviations from a lognormal density PDF can be 
estimated via the skewness Sg and kurtosis JCg of the logarithm of the density, s = ln{n/{n)), 
where (n) in the mean density. For a perfectly lognormal distribution, Sg = and JCg = 3. 
For the HD simulations, the deviations of the skewness and kurtosis from these fiducial values 



are h igher in the compressive forcing case than in the solenoidal forcing case (IFederrath et al. 



2OIOI . and Tabled]). We characterize the width of the density PDF by the density dispersion, 
cr„/(n), and the log-density dispersion, a^, which is the standard deviation of the logarithm 
of the density. The density PDF of the compressively forced field has a roughly three times 
higher standard deviation than the solenoidally forced counterpart at the same RMS Mach 
number, emphasizing the importance of studying different turbulent injection mechanisms. 
The values of o"„/(n) and a^ for the HD simulations are also given in Table [TJ 



3.2. Fractional Brownian Motion 



The hydrodynamic simulations have a unique velocity power spectrum, which cannot 
be varied. In other words, the compressively and solenoidally forced hydrodynamic simula- 
tions each provide one point in the calibration relation. In contrast, we aim to establish a 
relation between the PCA pseudo structure function and the slope of the energy spectrum 
over a range of exponents for the velocity spectrum, which reflects different types of turbu- 
lence (e.g., compressible, incompressible, intermittent, non- intermittent, etc..) in the ISM. 
This can be accomplished by varying the exponent of the energy spectrum, of the density 
spectrum, and t he dispersion of the density PDF of Fractional Brownian Motion structures 



(fBms, see, e.g., IStutzki et al.lll998l ). and by comparing the energy spectrum of the fBms 



to the PCA pseudo structure function obtained from the corresponding simulated spectral 
map. 



A comprehensive study of fBms by IStutzki et al.l (119981 ) summarizes the current knowl- 
edge of fBms, and we do not repeat it here. The details of the method we use to generate 



such fields are presented in lOssenkopf et al.l (120061 ) . and are summarized by the following. 
An fBm can be generated in Fourier space by creating an isotropic amplitude following a 
power law, A(k) = Aq k~'^. A phase 0(k) is randomly generated, using a uniform distribu- 
tion between — vr and tt to obtain the final Fourier transform of the desired field (velocity or 
density), /(k) = e^'^^^\ To ensure that the final field, obtained by taking the inverse Fourier 
transform of /(k), is real, the condition 0(k) = — 0(— k) is imposed. 

Twelve velocity fields and fifteen density fields were created on a 257^ grid with a power- 
law Fourier spectrum of exponents ranging between /3^ = 1.2 and /3^ = 3.4 for the velocity 
field, and /3„ = 0.6 to 3.4 for the density field. The exponents of the density and velocity 
fields were varied independently. This range of values for the velocity spectrum covers differ- 
ent types of turbulence, including Kolmogorov, Burgers, and log-Poisson tur bulence models. 
It also covers the case of systematic motions, such as infall, for which /3^ > 3 f lBrunt fc Heyer 
2002al ). Note that an amphtude Aq = 1 was used, and the final velocity and density fields 



were rescaled to a velocity standard deviation of 1 km s~^ and a mean density of 500 cm~^ a 
posteriori, as for the hydrodynamic simulations. For the density field, the rescaling is not as 
straightforward as for the velocity field, owing to its positivity. Hence, we subtracted from 
the original fBm its minimum value. The density field was then obtained by dividing the 
the fbm by its mean value and multiplying it by the desired mean density of 500 cm^^. The 
density fBm fields created with this method approximately follow a Gaussian (also known 
as normal) distribution, with a standard deviation of 90-180 cm~^. The density dispersions 
of each fBm density field created with a gaussian PDF and a varying density spectrum are 
listed in Table [2J 

Simulations and observations show that the density PDF of isothermal supersonic tur- 
bulent flows is better approximated by a lo gnormal distribution (i.e . , a Gaussian distri- 
bution in the logarithm of the density, see IVazquez-Semadenil Il994j: iPadoan et al.l 119971 : 



Passot fc Vazquez-Semadenilll998l : lFederrath et al.ll2008l : lPrice et al.ll201ll ). Furthermore, we 
wish to examine the dependence of the calibration relation not only on the exponent of the 
density Fourier spectrum, but also on the density PDF. In order to create d ensity fBms with 



a logn orma l PDF of variab l e stand ard deviation, the method presented in lOssenkopf et al. 



(120061 ) and iBrunt &: Heyerl (l2002al ) was followed. First, an fBm field was created with a 
power spectrum of slope /3„ = 1. This field represents the logarithm ln(?T,) of the desired den- 
sity field n. The fBm field ln(?T,) was then rescaled given the desired mean (n) and standard 
deviation o"„ of the desired density field. This rescaling is based on the relation between the 
mean and standard deviation of a lognormal field and its logarithm: 



(ln(n)) =ln((n))--af„(„). 



(1^ 
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0"ln(n) 



\ 




(2) 



where (ln(n)) and o'in{n) are the mean and standard deviation of the logarithm of the desired 
field. Thus, if (F) and ap are the mean and standard deviation of the fBm field F = ln(?T,), 
then F is rescaled and exponentiated to produce the final density field with the desired 
lognormal distribution: 



n = exp[{F- (F)) x ^i^ + (ln(n)) 



(3) 



Six density fields with lognormal distributions were created, with a mean of 500 cm ^ and 
standard deviations of 100, 1000, 2000, 3000, 4000, 5000, 6000, 7000, and 10000 cm'^ (a„/ < 
{n) ~ 0.2, 2, 4, 6, 8, 10, 12, 14, 20). Due to finite numerical resolution and low number 
statistics, the process of exponentiation may introduce small deviations from a lognormal 
PDF. To check the magnitude of the deviations of our fBms' density PDFs from a purely 
lognormal PDF, we computed the skewness and kurtosis of the logarithm of the density, 
which are listed in Table [31 The skewness is of the order of 0.01-0.04, so deviations from a 
purely lognormal distribution {Ss = 0) are much smaller than for the HD simulations [Sg = 
—0.1 and —0.26 from solenoidal and compressive forcing respectively). The kurtosis is also 
close (within 1.5%) to the value of /C^ = 3 obtained from pur ely lognormal distribu tions. 
Last, exponentiation changes the index of the Fourier spectrum ( jOssenkopf et al.ll2006[ ). The 
spectral indices of the density spectrum of the fBm created with a lognormal density PDF are 
listed in Table [31 The PDFs of the fBms generated via exponentiation are shown in Figure [H 
The dashed-lines represent the best lognormal fits to each PDF, the density dispersion of 
which is shown in the legend. 



3.3. Generation of the spectral maps 

A simulated spectral map of the ^^CO emission line (i.e., a PPV cube) was created for 
each simulation assuming that the ^ ^CO line is optically thin, and assurning a. n abundance 
ratio n(^3co)/n(H2) = 1.7 x lO"*' Jhaneer fc PenziasI Il990l : Isiake et al.lll987h . Note that 
the abundance ratio used in these simulations does not affect our results since it scales the 
CO intensity up or down but does not change the power spectrum of the density or velocity 
fluctuations. Thus, the simulated ^^CO PPV cubes were constructed from the density fleld 
n{x,y,z) and velocity flelds Vx{x,y,z), Vy{x,y,z), Vz{x,y,z) using the following expression 
(which is an example along the 2;-direction) along each cardinal direction [x, y, z): 
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/^(r) = ^j^(r,M)rf2; (4) 

z 

where r = (x, y, z), -u is the velocity channel, a^ is the absorption coefficient at the center 
of the ^'^CO line, By is the Planck function, and ly the specific intensity. The emissivity jy 
= ay By has units of W/m/str/Hz and is expressed by: 



/ N, -8 / 5.28 K 

jy{r,u)dz = 3.6x10 exp 



X 



T(r) 
K n(r) dz 



0.378 Tcm-3pcV2^zya,(r) 

e^^[-{u-v,{r)f/2a,{rY] 



where T is the temperature in K, n the number density field in cm~^, Vz is the projection 
along the z-axis of the velocity vector v at position {x,y,z) in ms~^, and the velocity 
dispersion a^ in ms~^, given by 

In this expression, "mco is the mass of the ^'^CO molecule. The first term represents thermal 
motions, and the second gas flows. Note that both the fBms and the hydrodynamic simula- 
tions are isothermal, with a uniform temperature of 10 K. In the following, we assume that 
the simulation box is L = 10 pc in size, and therefore, dz = L/N, where A^ is the number of 
grid points on one side (A^ = 256 for the HD simulations, A^ = 257 for the fBm simulations). 

The previous equations are based on two limiting assumptions: 1) the CO line is optically 
thin, and 2) local thermodynamic equilibrium (LTE). While the latter is true for densities > 
100 cm~^, CO is sub-thermally excited for smaller densities. Hence, the emission in the most 
diffuse regions of the simulations will be overestimated by our first assumption. On the other 
hand, CO becomes optically thick at column densities greater than N{CO) > 10^^ cm~^, 
which is not accounted for by our simple radiative transfer model. Such column densities 
can be attained in the high density regions of the simulations for a reasonable cloud depth 
(a few pc). However the scaling of the column density is arbitrarily set by the choice of the 
size of the simulation box, so it is pointless to try to determine whether this limit is actually 
reached in the simulations. 

The spectral maps resulting from hydrodynamic and fBm simulations were sampled on 
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a 40 m s~^ grid as in iFederrath et al.l (120 lOf ). To test the effects of spectral resolution on the 
uncertainty in the exponent of the PCA pseudo structure functions derived for each snapshot 
of the HD simulations or fBm field, we produced spectral maps with spectral resolution 10 m 
s~^ and 20 m s~^. Increasing the spectral resolution to 20 m s~^or 10 m s~^ did not reduce 
the scatter in the exponent of the PCA pseu do structure func t ion, a nd thus, we kept the 



original spectral resolution of 40 m s~^ used in IFederrath et al.l (J2010l ). A total of 486 ^^CO 
PPV cubes were generated from the hydrodynamic simulations (3 lines-of-sight directions 
for each of the 81 time snapshots for solenoidal and compressive forcing) and 255 ^^CO cubes 
were generated from the fBms (180 with variable density power spectra and 75 fBms with 
variable density PDFs). 



4. PCA applied to hydrodynamic simulations 

PCA was applied to all 486 PPV cubes generated from hydrodynamic simulations and 



a pow er law was fitted to each resulting PCA pseudo structure function. IFederrath et al. 



(I2OIOI ) presented the time-averaged PCA pseudo-structure function for the solenoidally and 
compressively forced hydrodynamic simulations, with slopes 0.66 ± 0.05 and 0.76 ± 0.09 
respectively. Here, we also derive the average slope of the PCA pseudo structure func- 
tion, averaged over all time snapshots and all three lines of sight (x, y, and z) and find 
(ttpcA)soi = 0.64 (standard deviation 0.05) and (apcA)comp = 0.77 (standard deviation 0.07) 
for the solenoidally and compressively forced simulations, respectively. Typical errors on 
the slope of the PCA pseudo structure function for individual snapshots are 0.02 and 0.04 
for solenoidal and compressive forcing respectively. In addition. Fig. |2] shows all the pairs 
of spatial and velocity scales detected by PCA in all the PPV cubes obtained from hydro- 
dynamic simulations. A power-law fit to this composite PCA pseudo structure function 
yields an exponent of apcA = 0.65 ± 0.05 and a'^^^ = 0.76 ± 0.07 for the solenoidally and 
compressively forced simulations, respectively. These results are in very good agreement 



with the slope of the time-averaged PCA pseudo structure function from IFederrath et al. 



( 120101 ). The results of PCA applied to individual time snapshots of the hydrodynamic simu- 



lations are shown as black crosses in Figure |3l their average being indicated by a red triangle. 
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5. PCA applied to fBms: Sensitivity of the PCA calibration to the density 

spectrum 

PCA was applied to the 180 PPV cubes generated from fBms velocity and density fields 
with varying power-law Fourier spectra. The colored lines in Figure [3] show the slope of the 
PCA pseudo structure function as a function of the exponent of the velocity spectrum for 
different density power spectra (the exponent of which is indicated in the legend). 

The slope of the PCA pseudo structure function, apcA, increase s with /3„, in agreemen t 



with previous calibrations based on fBms shown by the dashed line ( iBrunt fc Heyerll2002al ). 
The variation of apcA with /3„ is independent of the exponent of the density Fourier spectrum 
over the range /3„ = 0.6 to /3„ = 3.4. In addition, the average relation between ctpcA and 
/3„ obtained from HD simulations is too high (by a factor 2-3cr) compared to the relation 
obtained from fBms with gaussian PDFs. 

There are two major differences between the HD simulations and the fBms: 1) the 
HD simulations are intermittent, while the fBms are not; and 2) the fBms and HD density 
fields have different density PDFs, both in shape and standard devi ation. The fBms have a 



gauss ian density PDF of density dispersion anl{n) ~ 0.2-0.3, while iFederrath et al.l ( 12008 



20101) showed that the density PDFs of the HD simulations approximately follow lognormal 
distributions (i.e., gaussian in the logarithm of n). The discrepancy between the PCA pseudo 
structure functions obtained from HD simulations and fBms therefore suggests that the PCA 
calibration depends on the level of intermittency (both in the velocity and density fields), on 
the shape of the density PDF, and its density dispersion. This hypothesis is further explored 
in the next section. 



6. Sensitivity of the PCA calibration to intermittency and to the density PDF 

6.1. PCA applied to fBms with lognormal density PDFs of varying dispersion 

We have further tested the variation of the PCA calibration relation with density PDF 
by applying PCA to 180 PPV cubes generated from fBms with lognormal density PDFs 
of varying standard deviations, ranging from 100 to 10000 cm~^ (o"„,/(r7,) = 0.2-20, ag = 
0.2-2.45). Since the fBm velocity fields are not intermittent, we can thus isolate the effects 
of the density dispersion independent of the effects of intermittency in the velocity field, 
manifest in the HD simulations. The fBm density fields with varying density PDFs were cre- 
ated according to the method described in Section 13.21 Although a constant spectral index 
/Sn = 1 ch aracterizes the logarith m of the density fields, exponentiation changes the power 



spectrum (lOssenkopf et al.l 120061 ). As a result, the spectral indices of the fBm density fields 
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with lognormal PDFs are not equal to /3„ = 1, and are listed in Table [31 Nonetheless, we 
have shown that the calibration relation is insensitive to the index of the density spectrum 
in the previous section. Thus, the variations of Pn in the density fields with lognormal PDFs 
should not cause any variations in the calibration relation. 

The colored lines in Figure H] show the PCA calibration obtained for lognormal density 
PDFs of varying standard deviation (indicated in the legend). At high /3^, the PCA calibra- 
tion becomes unstable because it depends very strongly on a few Fourier components. We 
do not take into account values of (3^ > 2.6 in the following. For cr^ < 2, we do not find 
any significant variation in the PCA calibration as a function of log-density dispersion, while 
there is a sudden increase in apcA and its scatter for a given /3^ for a^ > 2. Actually, the 
PCA calibration becomes quite unstable for o"s > 2. For as < 2, the average calibration is 
shown by the black solid line. We derive a linear fit to the average PCA calibration obtained 
from fBms with lognormal density PDF of dispersion ag < 2, valid in the range /3^ = 1.2-2.6: 

(3^ = 0.20 ± 0.05 + (2.99 ± 0.09)apcA (6) 



This calibration is essentially identical to the relation derived in iBrunt fc Heyerl (J2002al ). 
within the errors. 

The increase and instability in apcA for cr^ > 2 is likely due to the inability of the 
density field to properly sample the velocity field at such high density dispersion. Extreme 
density fluctuations intermittently sample the velocity field, producing an effect similar to 
intermittency in the velocity field itself, similarly to discontinuous velocity jumps. In fact, 
we attempted to perform a PCA run on an fBm with a.^ = 3, but the field was so extreme 
that no scales could be detected in the PCA pseudo structure function. Note that, although 
the density dispersion of the compressively forced simulations ((T^°™p = 3.04) is higher than 
the density dispersion of the fBm with a^ = 2.45, this effect is not as strong for the HD fields 
because a high dispersion exponentiated fBm field is dominated by a few very high density 
point-like structures, while the hydrodynamically-produced density fields are dominated by 
a collection of filament-like structures. The latter are more spatially coherent, and capable 
of (~uniformly) sampling the velocity field across a longer region of space than the former. 
This effect, examined in more details in the next section, also appears to be a threshold 
effect, as shown by the absence of significant variations in the PCA calibration for o"s < 2. 



6.2. Effects of the density dispersion 

The PCA calibration derived from fBms with lognormal density PDFs of varying density 
dispersion exhibits some dependency to the log-density dispersion, cr^, but only above the 
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threshold o"s > 2 (see Fig. H]). The comparison between the PCA cahbration derived from 
fBms with cr^ < 2 and HD compressively forced simulations, which have a log-density dis- 
persion cr^°™p = 3.04 above the cr^ ~ 2 threshold, also supports the hypothesis that the PCA 
slope depends on cXs for a given /3„. Indeed, in Fig. HJ the PCA slope of the compressively 
forced simulations, shown as an open diamond, is apcA = 0.76±0.07, while the calibration 
relation obtained from fBms of density dispersion below the threshold cr^ < 2 predicts apcA 
= 0.58±0.03 for the corresponding Py = 1.94±0.05, or a factor of 3a difference. In order to 
prove that this difference in PCA slope is indeed due to the higher density dispersion of the 
HD compressively forced density field, we have performed two tests. 

First, we have generated spectral maps with the non-intermittent fBm velocity field of 
velocity spectrum /3„ = 2 and the 81 snapshots of the compressively forced HD density field. 
In this manner, we can isolate the effects of the HD density field from the effects of intermit- 
tency in the HD velocity field. The average PCA slope of all snapshots is shown as an open 
square in Fig. |H This PCA run with the fBm velocity field and the HD density field can 
reproduce, well within the errors, the PCA slope of the spectral maps generated from HD 
density and velocity fields with compressive forcing (open diamond), and demonstrates that 
the log-density dispersion of the HD density field contributes significantly to the increase in 
PCA slope compared to the average PCA calibration obtained from fBms with ag < 2. 

Second, we have rescaled both the HD compressively and solenoidally forced density 
fields to several log-density dispersions. For the compressively forced simulations, which 
originally have a density dispersion above the variation threshold of 0"^ — 2, we have scaled 
the log-density dispersion down to values (cr„/(n) = 0.8, 1.4, 2.0, 3.4, or 0"^ = 0.8, 1.3, 1.7, 
2.3) below or around this threshold in order to test if the PCA slope can be decreased down 
to values consistent with the PCA calibration obtained from lognormal fBms with ag < 2. 
For the solenoidally forced simulations, which originally have a density dispersion below the 
variation threshold of cr^ ~ 2, we have scaled the log-density dispersion up to values (cr„/(n) 
= 3.6, 5.9, 7.1, 11.8, or o"s = 1.9, 2.3, 2.4, 2.8) around or above this threshold. In this case, we 
thus test whether, according to expectations, the PCA slope obtained from the HD rescaled 
density fields with solenoidal forcing increases up to the level of the original compressively 
forced HD simulations for as above the variation threshold of the log-density dispersion. 

The rescaling was done similarly to the fBms with lognormal density PDF described in 
Equation [31 i.e. by rescaling the log of the density {F in Equation [3]) to the desired mean 
and dispersion calculated from Equations [1] and [2] and from the desired rescaled o"„/(n). 
We then exponentiated the rescaled log-density field. Because the HD density fields deviate 
from a lognormal distribution, the resulting log-density dispersion ag is not exactly related 
to the input o"„ by Equations [T] and [H but these equations nonetheless provide a good ap- 
proximation. The exact rescaled values of o"„ and o"s cited above are derived directly from 
the rescaled HD density fields. 
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Finally, we applied PCA to the spectral maps generated from the rescaled HD density 
fields and the non-intermittent fBm velocity field with /3^ = 1.9 for the solenoidally forced 
density field, and (3y = 2 for the rescaled, compressively forced HD density field. Again, the 
HD rescaled density fields are combined with non-intermittent fBm velocity fields in order 
to isolate the effects of the density dispersion from the effects of intermittency in the HD 
velocity field, which can potentially affect the calibration, as shown in Section 16.31 Fig. [5] 
shows the variations of the difference between the resulting PCA slope, apcA, and the PCA 
slope predicted from the calibration obtained with fBms of density dispersion a^ < 2, Op^^^ 
(see Equation [6]), as a function of as- In this Figure, we have also included the variations of 
C(pcA — C(pcA obtained from fBms with f3y = 1.8 and /S^, = 2.0. The shaded area indicates 
the 1, 2, 3 cr uncertainty in ap^^ from darkest to lightest. Note that apcA — ctpcA i^ not 
identically zero for fBms with cXs < 2 because the calibration derived in Equation E] is a fit to 
the average trend obtained from fBms with a^ < 2. Fig. E] demonstrates that, for ag < 2, the 
calibration is constant with cTs within the errors, while apcA starts to increase and deviate 
significantly from the calibration for cxs > 2 due to the poor sampling of the velocity field 
by the density field. Above o"s ~ 2, the variations of apcA with as for a given (3^ are steeper 
and more uncertain for the fBms than for the HD rescaled fields, which is also seen in Fig. 
m As mentioned before, a probable explanation for this difference is that the fBms are just 
dominated by a few point-like structures that sample the velocity field very poorly, while the 
HD density field structure consists of filament-like features which are more spatially coherent 
than the fBm density extrema. 



6.3. Effects of intermittency in the velocity field and the operating order of 

PCA 

Although the log-density dispersion of the solenoidal forced simulations (cr^°' = 1.3) 
is below the threshold of cTs ~ 2 above which apcA starts to increase with as, Figure H] 
shows that the PCA slope of the solenoidally forced HD simulations still stands out as being 
too high compared to the calibration derived from fBms. Indeed, the average slope of the 
PCA pseudo structure function is 0.65±0.05, while the average calibration obtained from 
fBms with density dispersion cxs < 2 predicts apcA = 0.55±0.03 for the corresponding /3„ = 
1.86±0.05. In order to determine the cause of this discrepancy, we have first confirmed that 
the HD solenoidally forced density field was not causing this difference by applying PCA to 
spectral maps generated from a non-intermittent fBm velocity field of velocity spectrum P^ 
= 1.9 and the 81 snapshots of the solenoidally forced HD density field. The result is shown 
by the open triangle in Fig. H] (this point is also shown in Fig. [5]). The average PCA slope of 
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the spectral maps generated from the solenoidally forced HD density fields and fBm velocity 
field with /3j, = 1.9 is in nearly perfect agreement with the PCA calibration derived from 
fBms velocity and density fields with lognormal density PDF of dispersion 0"^ < 2. In Fig. El 
the combination of the HD solenoidally forced density field and the fBm velocity field with 
I3y = 1.9 fits well within the errors in the trend apcA — Ci'pcA versus ag. This demonstrates 



that the HD solenoidally forced density field is not causing the PCA slope obtained for HD 
solenoidally forced simulations to be too high compared to the calibration obtained for the 
same range of log-density dispersion. We conclude that, in this case, this discrepancy must 
be due to the intermitt ent structure of the HD velocity field. 



Brunt et al.l ( 120031 ) showed that the slope of the PCA pseudo structure function reflects 
the slope Ci of the first order structure function, defined as SFi(£) = (t>(r + i) — v{r)), and 
that the calibration between ({1) and apcA is insensitive to the level of intermittency. On 
the other hand, they showed that the relation between apcA and the slope C(2) of the second 
order structure function, defined as SF2(£) = {[v{r + i) — v{r)] ), or equivalently the slope of 
the velocity spectrum / 3„ = 2C(2) + 1, d epend s on the level of int ermittency. It was de mon- 
strated theoretically by lBoldyrev et al.l ( 12002| ) and numerically by lSchmidt et al.l ( l2008l ) that 
the relation between the order p of the structure function SFp{£) = {[v{r + i) — v{r)Y) and 
its exponent ({p) is more concave as the level of intermittency in the velocity field increases 
(i.e., it rises slower than linear). Thus, the ratio of C(l)/C(2) increases as velocity fields 
becomes more intermittent. Because PCA traces C(l)) it is therefore expected that apcA 
increases with intermittency in the velocity field for a given /3„. This is confirmed by the fact 
that the intermittent HD solenoidally forced velocity field (rather than the density field) is 
causing the increase in apcA compared to the fBms of same (3y (and log-density dispersion), 
which are not intermittent. In the following, we quantify this effect for the HD simulations. 

The level of intermittency in the velocity field increases from fBms (non-int ermittent), to 
soleno idally forced, to compressively forced HD simulations, as demonstrated in lSchmidt et al. 
( I2OO8I ). As a result, the ratio C(l)/C(2) is lowest for the fBms (which are not intermittent 
and for which ({p) is linear with p, see iBrunt et al.ll2003l ). increases for solenoidally forced 
HD simulations, and increases even more for co mpressively for c ed sim ulations. We can esti- 
mate the ratio C(l)/C(2) using the results from lSchmidt et al.l ( 120081 ) . who found ({l)/({3) 
= 0.47±0.01 and C(2)/C(3) = 0.79±0.01 for the case of solenoidal forcing, and C(l)/C(3) = 
0.63±0.01 and ({2)/({3) = 0.90±0.01 for the compressively forced HD simulations. We can 
then derive ({1) for each forcing case using the following Equation: 



C(i) 



rc(i)i 

[C(3)J 


rc(3)i 
Lc(2)J 



(/3. - 1) 



(7) 



where we used the fact that C(2) = /3y — 1. We find ({1) = 0.51±0.03 for the solenoidal 
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forcing case and_C(l) = 0.66 ±0.04 for the compressive forcing case. Using Equation 31 in 



Brunt et al.l ( 120031 ) . which relates ({1) and apcA, with 7 = ^(1) (the operating order of 
PCA), we predict apcA = 0.62±0.03 for solenoidal forcing and apcA = 0.75±0.07 for com- 
pr essive forcing . This prediction, based on the relation between ^(1) and apcA established 



m 



Brunt et al.l (120031 ) and the characterization of intermittency by the ({p) relation in the 
HD velocity field, is consistent within the error bars with the PCA slopes derived here for 
the HD simulations {af^A = 0.65±0.05 and a™^^ = 0.76±0.07). 

The difference between apcA obtained for HD simulations and fBms of similar log- 
density dispersion and velocity spectrum (see Fig. H]) therefore confirms that intermittency 
in the velocity field changes the relation between I3y and apcA- For a given (3y, the variations 
of a PC A between HD and fBms simulations can then be simply explained in terms of 1) the 
operating order of PCA and 2) the variations of the relation ((p) between the exponents of 
the structure functions of different orders p with intermittency in the velocity field. In fact, 
one goal of this paper is precisely to investigate how the calibration between apcA and /3„ 
varies with intermittency and how to account for it in the derivation of (3^ from PCA. For 
the solenoidal forcing case, the difference in PCA slope between non-intermittent fBm and 
intermittent HD fields is entirely due to intermittency in the HD velocity field, since, in this 
case, the log-density dispersion of the density field is below the as — 2 threshold for which 
poor sampling of the velocity field causes the PCA slope to increase at constant /3„. For the 
compressive forcing case, we have demonstrated that the lack of sampling of the velocity 
field by the density field due to the high log-density dispersion contributes significantly to 
the difference in PCA slope between the compressively forced HD simulations and the fBms 
of same velocity spectrum. Nonetheless, the contribution from intermittency in the velocity 
field likely plays an important role, although it is difficult to detect it considering the larger 
error bars obtained for the HD simulations with compressive forcing. 



6.4. Consequence for molecular cloud observations 

Applying a PCA calibration obtained from the idealized structures that are fBms to 
actual molecular clouds with high levels of intermittency and very high density dispersions 
will lead to an overestimation of (3^. For instance, in the case of the solenoidally forced 
simulations {apcA = 0.65±0.05, 5, = -O.lOiO.ll, /C^ = 3.03±0.17), the PCA calibration 
obtained from fBms with cr^ < 2 predicts (3^ = 2.12±0.17, while the actual value of /3„ 
measured in the simulated velocity field is /3^, = 1.86±0.05. For the compressively forced 
simulations {apcA = 0.76±0.07, 5, = -0.26±0.20, /C^ = 2.91±0.43), the PCA calibration 
obtained from fBms estimates /3„ = 2.48±0.27, while actually, /3„ = 1.94±0.05. The difference 
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and induced uncertainty in these values of (3^ appears very large compared to the subtle 
difference between, e.g., Kolmogorov turbulence {(3^ = 5/3) and Burgers turbulence {(3y = 
2). As a result, one should explore values of /3^ implied by molecular line observations 
based on PCA calibration relations derived for different density dispersions and levels of 
intermittency in the velocity field. 



7. Application to molecular clouds identified in the Galactic Ring Survey 

In this section, we apply PCA to spectral maps of 750 molecular clouds identified in 
the Five College Ra dio Astronomy Observatory (FCRAO) Galactic Ring Survey (GR S, see 



Jackson et al. 



20061 ) for which distances are available from iRoman-Duval et al.l (120091 ). Sev- 



eral observational difficulties must first be overcome before PCA can be applied to observa- 
tions of molecular clouds. 



7.1. Subtraction of the noise ACF 



The noise inherent to observations of molecular clouds contributes to the observed ACF 
of the principal components. It can be shown that the no ise in the principal components 
is identical to the noise in the data (JBrunt fc Heyerll2002al ). For spectroscopic data cubes 
obtained from on-the-fly mapping (as in the GRS), the noise is correlated between positions 
observed with the same reference position. As a result, the ACF of the noise contributes a 
powerful peak at the z ero-lag, and other pe aks corresponding to different correlation lengths 
of the noise (see also I Jackson et al.l 120061 ). The ACF of the noise therefore needs to be 
subtracted from the observed ACF of the principal components. The contribution of the 
noise to the ACF of the principal components was estimated by averaging the ACF of the 
5 highest-order principal components, which should only contain contributions from noise, 
because they correspond to features of spatial scales well below the resolution limit. The 
ACF of the noise was then subtracted from the ACFs of all the principal components before 
computing the spatial scales. 



7.2. Effects of the finite telescope resolution on the determination of spatial 

scales 



The convolution of the observed ^^CO data with the beam of the telescope affects the 
determination of spatial scales. The removal of the beam contribution to the observed spatial 
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scales has been established by iBruntI ( ll999l ): lBrunt fc Heyerl (l2002al ). using approximations. 
Here, we analytically derive the exact scale correction for beam smearing, and test it using the 
HD numerical simulations. Let X be a two-dimensional image observed through a telescope 
with a two-dimensional beam B, and let X be the "true" image (i.e., not convolved with the 
telescope beam). X is the convolut ion of the tr ue image with the telescope beam: X = X-kB. 
It can be demonstrated (see, e.g., lBruntlll999l ) that the ACF of the beam-smeared principal 



components, ACF — ^(„), is related to the ACF of the true principal components, ACFp(-,(n), 



by ACF^(„) 



ACFpc(n) * ACFb, where ACF^ is the ACF of the beam. The FCRAO 
telescope beam is well approximated by a Gaussian beam of width gb'- 



B(r) 



27ra| 



exp 



2al 



such that the corresponding normalized ACF is 



ACF B{i) = exp 



4a| 



(9) 



The ACF of a turb ulent field and its principal compon ents can be well approximated 
by an exponential ACF jBruntlll999l : iBrunt fc Hevei1l2002ah : ACFpc(n)(£) = e"^, where 
A*^"^ is the true spatial scale corresponding the n*'^ principal component. This is shown in 
Figure [6l where the crosses represent the (true) ACF of the principal components of one time 
snapshot of the hydrodynamic simulation (with solenoidal forcing). The solid line indicate 
the exponential function with the same e-fold length as calculated in the ACF (indicated in 
units of pixels in the legend). An exponential function fits the ACF well, and seems to be a 
reasonable assumption. Note that only the core of the ACF (above the 1/e level) matters for 
the calculation of the e-fold length, and thus the fact that an exponential function does not 
fit the wings of the ACF is unimportant. Combining the exponential form of the true ACF 
of the principal components and the relation between the ACF of the observed principal 
components and the true principal components, the un-normalized ACF of the observed 
principal components is: 



ACF™,„,(£) 



exp 



x=0 



Xin) 



exp 



(^ - X)2 . , 
2 ) "^+ / 6^P 



4a 



x=0 



X 
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+ x 



4a| 



dx 



(10) 



The two integrands, /_ and /+, can be factored as: 
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I± = exp 



4a| 



x + 2{ 
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exp 
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We then make the change of variable y = (x + 2(o-|/A'^"'^) ± ^/2)) /(2aB), and we recognize 
the error function (erf): 



e-^'rfi/ = ^(l-erf(a)) 



(12) 



so we get: 



ACF;1(„)(£) 



y/naB < ( 1 — erf 
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(13) 



Then we estimate ACF'l!l(„)(0): 

PC 



ACF™ („, (0) = 2V^aB (^1 - erf (^ ) ) exp 



o-B x2 



L^AW- 



(14) 



and finally get the ACF of the observed principal components normalized so that ACF~(„) (0) 
1: 



ACF~(„)(£) = —, —^ 

P^ ^^ 2(l-erf(^)) 



1 — erf 



1 — erf 
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(15) 



In order to relate the observed and true spatial scales, we need to determine for what value 
of i the auto-correlation function of the principal components, ACF~(„)(£), falls by one 
e-fold as a function of the true spatial scale \^^\ For true spatial scales A*^"^ ranging from 
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5" to 5000" (sampled every 5"), we computed ACF~,(n)(£) according to Equation [T5] and 

determined the corresponding "observed" e-fold length X^^l. We then constructed a look-up 
table relating A*-"-* and X^^i., which is the desired scale correction. For every spatial scale 
detected in the GRS molecular clouds, the scale correction was applied by finding the two 
closest X^^ values in the look-up table, and interpolating the corresponding true scale A^"^ 
accordingly. 

We used the hydrodynamic simulations presented in Section 13.11 to test our scale cor- 
rection. We convolved one of the simulated ^^CO data cubes generated from hydrodynamic 
simulations with a Gaussian beam of FWHM width 48", and applied PGA to the resulting 
beam-smeared simulated PPV cube, with and without scale correction. The result is shown 
in Figure [71 The black crosses represent the scales detected in the simulations without convo- 
lution with the beam (i.e., the true scales). The blue triangles represent the scales detected 
in the simulations after convolution with the beam, but without any scale correction. Fi- 
nally, the red stars indicate the scales detected in the simulations convolved with the beam 
after scale correction. For comparis on, the green diamonds show the scales corrected with 



the (approximate) prescription from lBrunt &: Heyerl (J2002al ). The power- law nature and the 



exponent of the PGA pseudo structure function is only recovered after correcting the spatial 
scales for the convolution with the telescope beam. The exact analytical scale correction 



present ed here provides more accurate results than the prescription from iBrunt fc Heyer 



( l2002al ). which is based on an approximation. Last, we emphasize that, when applying this 
scale correction to real observations of molecular clouds, only spatial scales above the res- 
olution limit (as detected before the scale correction) must be taken into account. Spatial 
scales below the resolution limit correspond to noise and must be excluded from the analysis. 
In the following, we choose to use the FWHM of the beam as the resolution limit. This is 
justified in Section [m The changes in PGA slopes caused by different choices of resolution 
limit (e.g., 3cr, 2 x FWHM) are also investigated in Section E3 



7.3. Uncertainties in the PCA results 

The uncertainty in the spatial scales detected by PGA stems from the uncertainty 
in the distance. The error of the kinematic distances of the GRS clouds were estimated 
in 



Roman-Duval et al.l ( 120091 ) and are propagated here. The finite size of the pixels also 
contributes to the uncertainty on the spatial scales. Specifically, the error on the spatial 
scale is given by aj = aj x 69^ + {6pij./2 x d) , where ad is the error on the distance, Opi^ is 
the angular size of a pixel, and 59 is the angular scale detected by PGA. The uncertainty on 
the velocity scales was set to half the velocity resolution (0.1 km s~^ for the GRS). 
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7.4. Results 



PCA was applied to molecula r clou ds identified in the GRS for which distances are 



available from iRoman-Duval et al.l (|2009[ ). Out of the 750 molecular clouds for which dis- 
tances were available, 383 did not exhibit a large-enough spatial dynamic range to allow 
the detection of five or more spatial and velocity scales. Consequently, a robust power law 
could not be fitted to the resulting PCA pseudo structure function for this sample of clouds. 
Our sample of GRS clouds thus contains 367 objects. Figures [8] and [9] show an example of 
PCA results for a particular cloud, GRSMC G053.59+00.04. In Figure El the O*'^ principal 
component simply shows the integrated intensity of the cloud and thus provides information 
on its overall structure. In Figure IHl only spatial scales above the resolution hmit (before 
scale correction) are shown. A power law was fitted to the PCA pseudo structure function, 
yielding apcA = 0.74±0.05. 

The black histogram in Figure [10] shows the histogram of the slope of the PCA pseudo 
structure function obtained from the GRS clouds. The mean PCA slope is (apcA) = 
0.61 ± 0.2, where the error bar reflects the standard deviation of the distribution. To reduce 
the effects of outliers, we also computed the average of apcA weighted by the inverse of the 
reduced x^ of the power-law fit to the PCA structure function, and obtained a weighted 
average {apcA)w = 0.62±0.2. This value of apcA is in good agreement with the PC A slope 



obtained for molecular clouds located in the Outer Galaxy ( iBrunt fc Heyerl l2002bl . (apcA) 
= 0.62±0.11). A power law of slope 0.62 also fits well the composite structure function, 
composed of all the spatial and velocity scales detected in all the clouds (see Figure [TT]) . 
A bisector fit to the PCA composite structure function shown in Figure [11] yields apcA 
=0.60±0.2. 



7.5. Estimation of the density dispersion of GRS molecular clouds 



Since the measured value of apcA is unstable above as ~ 2, it is worthwhile trying 
to estimate plausible values of as that may be present in the GRS cloud sample, to gauge 
the possible effects of high density dispersion on our results. Models of driven turbulence 
suggest that the density dispersion is related to the 3-dimensional rms Mach number (A^) 
as follows: 0"„/( n ) = bAi, where 6 is a c onstant dependin g on the nature of the turbulent 
driving . For solenoidal forcin g, b ^ 1/3 (iPrice et al.ll201ll ). while for compressive forcing, 
6^1 (JFederrath et al.ll2008l ). There are very few obs ervationally-determined values of b, 
but existing measurements favor b ^ 0.5 (JBruntll2010l ). indicating a mixture of solenoidal 
and compressive forcing. 
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Assuming a lognormal PDF, so that as = \/ln(l + [an/ (n)y, and with a specified 
kinetic temperature, T and mean molecular mass, m, we can derive a relation between the 
1-dimensional velocity standard deviation, ct^^id, and as as follows: 
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(16) 



Here, we have used M. = Cu^sd/cs, where Cs = ^^/kT/m is the sound sp eed, and have assu med 
a mean molecu lar mass of 2 .72 times the ma ss of a hydgrogen atom (|Hildebrandlll983l ) and 
taken b = 0.5 (JBruntll2010l ) and T = 10 K (JRoman-Duval et al.ll2010l ) as reference points. 



The c hoice of a kinetic temperature of 10 K is motivated by Fig. 6 in iRoman-Duval et al. 



(J2010l ). where the maximum excitation temperature in a molecular cloud occurs in the dens- 
est regions that are closest to LTE, and should refiect the actual kinetic temperature of the 
gas. We have also assumed isotropy, so that the 3-dimensional velocity standard deviation 



V alues oia^^iD for the GRS cloud sample have already been measured by lRoman-Duval et al 
( 120101 ). We have converted these measurements into estimates of as for the sample of 367 



clouds analyzed here, and the histogram of the resulting as values is shown in Figure [T2l 
The histogram peaks near as = 2.1, with a tail extending to as ~ 2.4. Comparison of the as 
histogram with the HD results in Figure [5] suggests that a minor overestimation of apcA may 
be present in some clouds due to extreme density fiuctuations (up to ~ +0.1). In general 
though, as long as 6 = 0.5 and T = 10 K reasonably represent the conditions in the GRS 
clouds, then we conclude that extreme density fiuctuations have a relatively minor impact on 
our measured apcA- While it is unlikely that kinetic temperatures are below 10 K, if extreme 
compressive forcing is common then the as values will be a little higher than represented in 
Figure [12] (but note that as varies only slowly with b due to the square-root of a logarithm 
dependence and this may be countered by raised kinetic temperatures). 



7.6. Turbulent spectrum of GRS molecular clouds from PCA 

Due to the large dispersion of apcA, PCA provides a coarse measurement of (3y for any 
individual cloud. However, when considering the ensemble average, it is a reliable statistical 
measure of the exponent of the turbulent spectrum. Applying the PCA calibration derived 
in Section m based on non-intermittent fBms with as < 2, the mean value of the PCA slope 
(0.62±0.2) corresponds to (/3„) = 2.06 ± 0.6, where the error bar refiects the standard de- 
viation (see Figure [TOj) . The large standard deviation refiects not only the uncertainty on 
the derivation of apcA, but potentially also intrinsic variations of the turbulent spectrum 
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between different molecular clouds, due to varying star forming activities, different sources of 
forcing (e.g., solenoidal versus compressive) and driving scales, and a range of Mach ii umbers 
flKlessenll200l[ lOssenkopf k Mac L0JI2OO2I : ISrunt et allbood : IPederrath et al.ll201oh . 

As pointed out before, intermittency and high density dispersion can introduce signifi- 
cant deviations compared to predictions from fBms. Since these effects are likely to play a 
significant role in molecular clouds, as shown by the HD simulations, we need to also com- 
pare the results of PCA applied to GRS molecular clouds with HD simulations. The average 
slope of the GRS PCA pseudo structure functions (0.62±0.2) is in excellent agreement (< la) 
with the PCA slope derived from the spectral maps generated from solenoidally forced HD 
simulations {{apcA) = 0.65±0.05), and in marginal agreement (2cr) with the compressively 
forced simulations {{apcA) = 0.76±0.07), which likely exhibit a higher density dispersion 
than the average GRS molecular cloud, as demonstrated in Section 17.51 For these two cases 
of turbulence forcing, the exponents of the energy spectrum measured in the velocity fields 
are (i^ = 1.86±0.05 and /3„ = 1.94±0.05, respectively. Therefore, accounting for intermit- 
tency and density dispersion effects yields a spectral energy index of /3„ ~ 1.9 for the GRS 
molecular clouds. This range of value s correspond to log-Poisson (intermit t ent, compress- 



ible) turbulence (IShe fc Levequd ll994J : iBoldyrevI l2002l : iBoldyrev et al.l |2002[ ISchmidt et al. 
20081 ). but is also consistent with Burge rs turbulence w i thin t he errors. 



For the HD simulations used here, ISchmidt et al.l (|2008[ ) showed that the relation be- 
tween the scaling exponents ({p) of the structure functions of orders p = 1-5 are consistent 
with a log-Poisson model, for which C(p)/C(3) = (1 — '^)3 + C{1 — (1 — ^)^), where A and 
C are the scaling exponent (or second order structure function exponent) and co-dimension 
of the dominant dissipative structures, respectively. Both C and A depend on the degree of 
intermittency of the flow. For the HD simulations discussed here, and assuming A = 1 (the 
dominant dissipative structures are shocks that obey Burgers turbulence scaling relations), 
Schmidt et al.l (120081 ) found C = 1.1 and 1.5 for the compressive and solenoidal forcing cases 
respectively (i.e., the dominant dissipative structures are 2D shocks). In contrast. Burgers 
turbulence predicts an exponent C(p) = 1 for p >1, inconsistent with the scaling exponents 
of the structure functions in the HD simulations. Hence, if the HD simulations are an accu- 
rate model of molecular clouds, the coincidence between the exponents of the PCA pseudo 
structure functions derived from the HD simulations and the GRS molecular cloud catalog 
suggests that turbulence in molecular clouds is best described by a hierarchical, intermittent 
log-Poisson turbulence model with 2D shocks as the most dissipative structures. However, 
we cannot formally distinguish between log-Poisson models with 2D shocks as singular dis- 
sipative structures and Burgers turbulence in the GRS molecular cloud sample due to 1) the 
large uncertainty in the average exponent of the energy spectrum derived from the GRS cloud 
sample and 2) the fact that the velocity field and structure functions are not observable. 
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7.7. The choice of "resolution limit" 

Throughout Section [3, we applied a cutoff to spatial scales detected by PCA. Only 
spatial scales above the "resolution limit" as calculated before scale correction, with the 
resolution limit being defined as the FWHM of the beam (48"), were considered in the PCA 
pseudo structure functions of GRS molecular clouds. Scales below this threshold are consid- 
ered to be the result of noise and are excluded from the analysis. However, scales between 
the Nyquist sampling scale and twice the FWHM of the beam probably contain contribu- 
tions from both astrophysical signal and noise, and our decision to use the FWHM therefore 
needs to be justified. 

Figure [12] shows the composite PCA structure function, with the spatial scales left as 
angular scales in units of arcsecs uncorrected for beam smearing. Different characteristic 
values of the gaussian beam (la, 3cr, FWHM, 2 x FWHM) are indicated by vertical lines. 
The progressive loss of information between spatial scales corresponding to Nyquist sampling 
(which is equal to the la width of the beam, or 20") and the 3a width of the beam (60") 
is seen in the angular composite PCA structure function as a progressive change of slope, 
the slope becoming shallower as spatial scales get closer to the Nyquist sampling limit. For 
scales below the Nyquist sampling limit, the PCA pseudo structure function only reflects 
contribution from the noise, all velocity scales are equal to the spectral sampling (0.2 km 
s~^), and the PCA pseudo structure function is fiat (slope zero). This progressive decrease 
in PCA slope with decreasing spatial scales starts between 48" and 60" (the FHWM and 3a 
width of the beam). At spatial scales corresponding to 20" (the Nyquist sampling scale and 
the la width of the beam), the PCA angular composite structure function is dominated by 
the noise. Therefore, choosing the la width of the beam as the resolution limit would result 
in a mean PCA slope significantly skewed by the contribution of the noise. Hence, we took 
a conservative approach and excluded scales above the FWHM of the beam. 

Nonetheless, we have investigated the changes in PCA slope incurred by different defi- 
nitions of the resolution limit. The purple, red, and blue histograms in Figure [TO] show the 
distributions of apcA and /3y obtained from resolution limits of 20" (la beam width), 60" 
(3a beam width), and 96" (2 x FHWM) respectively. The resulting mean values of apcA and 
Pv become steeper as the resolution limit increases (going from la of the beam to its FHWM 
to 3a to 2 X FWHM). This is due to the decreasing contribution of noise as the definition of 
the resolution limit becomes more conservative. As expected, the mean PCA slope obtained 
from the la definition of resolution limit is significantly lower than for the 3 other cases 
because it includes a large contribution from the noise. Excluding this case, the mean values 
of apcA and (3^ calculated with different definitions of the resolution limit are within the 
errors of each other, and all consistent well within the errors with the HD simulations. We 
conclude that, as long as the resolution limit is above the FWHM of the beam, the choice of 
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resolution limit (between FWHM, Sa, and 2x FWHM) does not change the interpretations 
and the conclusions presented here. 



8. Summary and conclusion 

We applied Principal Component Analysis (PCA) to synthetic Position- Position- Velocity 
(PPV) spectral maps generated from the density and velocity fields of solenoidally and com- 



press i yely f orced hydrodynamic simulations of supersonic turbulence (JFederrath et al.l 12008 



20091 . I2OIOI ). and of fractional Brownian motion simulations, in order to constrain the cali- 
bration relation between the PCA pseudo structure function and the index of the velocity 
spectrum of turbulence, and to examine the dependency of this relation on the density spec- 
trum, intermittency, and density dispersion. 

We demonstrated that the calibration relation, the relation between the slope of the 
PCA structure function apcA and (3^, does not depend on the exponent of the power- law 
density spectrum /3„. 

For a log-density dispersion cr^ < 2, we do not find any dependence of the PCA 
calibration on the dispersion of the density PDF. We derive a PCA calibration relation, 
/3„ = 0.20 ±0.05 + (2.99 ±0.09)apcA valid for a, < 2 and /3^ = 1.2-2.6. For a, > 2, we find a 
strong dependence of the calibration between apcA and /3„ with ag- Extreme density fiuctua- 
tions intermittently sample the velocity field, producing an effect similar to intermittency in 
the velocity field itself - i.e. mimicking discontinuous velocity jumps, although the detailed 
mechanism is rather different. PCA is stable below a threshold of the log-density dispersion, 
as ~ 2, but if real molecular clouds exceed this, then an additional overestimation factor 
applies to apcA- Without knowledge of the true 3D log-density dispersions in the cloud 
sample, the estima tion of the turbulent spectrum in molecular clouds remains uncertain. 
Brunt et al.l (l2010bl Jal) developed a method to estimate the density PDF of molecular clouds 
based on the 2D power spectrum, the variance, and the PDF of the 2D column density, 
from which the 3D density PDF can be reconstructed, even in cases where the density PDF 
is not lognormal. However, this method requires high fidelity measures of column density 
such as extinction derived from 2MASS photometry of background stars and high spatial 
dynamic range. Therefore, it is not readily applicable to our set of data from the Galactic 
Ring Survey for which the spatial dynamic range for most clouds is limited. In addition, 
numerica l simulations pred ict a relation between the log-density dispersion and the Mach 
number ( jPrice et al.l 120111 ). but t his relation also depen ds on the relative contribution of 
solenoidal and compressive modes (JFederrath et al.ll2010f). An initial te st of the log- density 
dispersion - Mach number relation has been made ( Brunt et al.ll2010bl Jal). and this suggests 
that both solenoidal and compressive forcing are important and that density dispersions are 
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likely to be high enough that their effect on PCA is not insignificant. 

We demonstrated that intermittency in the velocity field also increases the PCA slope 
for a given velocity spectrum. This effect is due to a combination of the operating order of 
PCA (PCA traces the first order structure function exponent, C(l))) ^^^ the variation of 
the ratio between ({1) and /3„ with intermittency. Thus, if a first-order scheme is used to 
measure the second order exponent /3„, then some knowledge of the level of intermittency 
is required. By accounting for the level of intermittency, we were able to reconcile PCA 
measurements between non-intermittent fBms and the HD fields. 

We applied PCA to ^^CO spectr al maps of 367 molecular clouds identified in the Galac- 



tic Ring Survey (I Jackson et al.ll2006l ). We found that the average slope of the PCA pseudo 
structure function and the slope of the composite structure function, made of all the spatial 
and velocity scales derived in all the GRS clouds, are consistent with apcA = 0.62 ±0.2. Ap- 
plying the PCA calibration obtained from fBms with a^ < 2, the PCA slope obtained for GRS 
molecular clouds corresponds to an average turbulence spectral index of (/9) = 2.06 ± 0.6. 
However, we have shown that intermittency and density dispersion need to be taken into 
account. The average PCA slope obtained for the GRS clouds is in very good agreement 
with the PCA slope obtained from both solenoidally and compressively forced HD simula- 
tions, albeit in better agreement (at < 1 a) with the solenoidally forced HD simulations. 
This agreement suggests that turbulence in molecular clouds, as in the HD simulations, obey 
log-Poisson scaling relations (intermittent, compressible turbulence) with 2D shocks as the 
dominant dissipative structures. 
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Table 1: Turbulence statistics obtained for hydrodynamic simulations 



Symbol 


Description 


Solenoidal Forcing 


Compressive Forcing 


Pn 


Exponent of the density spectrum 


0.78 ±0.06 


1.44 ±0.23 


Pv 


Exponent of the velocity spectrum 


1.86 ±0.05 


1.94 ±0.05 


o-„/(n) 


Standard Deviation of n 


1.89±0.09 


5.86±0.96 


{s) 


Mean of s = ln{n/{n)) 


-0.83±0.05 


-3.40±0.43 


O-s 


Standard Deviation of s 


1.32±0.06 


3.04±0.24 


Ss 


Skewness of s 


-0.1±0.11 


-0.26±0.20 


ICs 


Kurtosis of s 


3.03±0.17 


2.91±0.43 



Table 2: Density dispersion of the fBm density fields with varying /3„ and gaussian density 
PDF 



(3n 


an/{n) 


0"s 


0.600000 


0.20 


0.21 


0.800000 


0.18 


0.19 


1.00000 


0.19 


0.21 


1.20000 


0.21 


0.22 


1.40000 


0.22 


0.23 


1.60000 


0.23 


0.25 


1.80000 


0.23 


0.26 


2.00000 


0.25 


0.28 


2.20000 


0.28 


0.32 


2.40000 


0.27 


0.33 


2.60000 


0.36 


0.39 


2.80000 


0.34 


0.40 


3.00000 


0.34 


0.40 


3.20000 


0.32 


0.41 


3.40000 


0.30 


0.41 
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<J./<n> = 0-2000 

CTy<n> = 2.000 

CT„/<n> = 4.000 

o„/<n> = 6.000 

aj<n> = 8.000 

CT„/<n> = 10.00 

oJ<r\> = 12.00 

<r„/<n> = 14.00 . 

CT./<n> = 20.00 



-5 5 

s = ln{n/<n>) 



Fig. 1. — PDFs, assumed to be lognormal, of the fBm density fields generated from ex- 
ponentiation. The dashed hnes represent the best lognormal fit to each PDF, the density 
dispersion of which is shown in the legend. 



Table 3: Moments of the density PDF and spectral index, /3„, for the fBms with lognormal 
density PDF and varying density dispersion, cr„ 



cr„/(n) (s) 



cr.« 



5. 



/C., 



/3n 



0.2 -0.02 0.20 0.03 3.02 1.04 

2.00 -0.81 1.27 0.04 2.95 0.63 
4.09 -1.42 1.68 0.01 2.98 0.37 
5.92 -1.83 1.90 0.04 2.97 0.22 

8.01 -2.11 2.04 0.02 2.97 0.16 
10.02 -2.33 2.15 0.04 2.97 0.22 
12.55 -2.47 2.23 -0.01 2.99 0.11 
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Fig. 2. — PCA composite structure function for the hydrodynamic simulations, containing 
all spatial and velocity scales detected by PCA in all PPV cubes corresponding to solenoidal 
forcing (top) and compressive forcing (bottom). The solid line indicates the best fit. In 
the top panel (solenoidal forcing), the dashed line shows the best-fit from the compressive 
forcing case. In the bottom panel (compressive forcing), the dashed line shows the best-fit 
from the solenoidal forcing case. 
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Fig. 3. — Calibration relation between the slope of the energy spectrum, Py, and the slope 
of the PCA pseudo structure function, apcA, derived from PCA applied to fBms and hydro- 
dynamic simulations. The top and bottom panels correspond to solenoidal and compressive 
forcing, respectively. The colored lines correspond to the relation between /3j, and apcA ob- 
tained from fBms with different density spectra, indicated in the legend. The black crosses 
correspond to the calibration relation deduced from each hydrodynamic time snapshot, the 
average of which is shown by the red triangle. 
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Fig. 4. — Calibration relation obtained from hydrodynamical simulations and density fBm 
fields with lognormal distributions. The black star and triangle represent the average result of 
PCA applied to each time snapshot of the solenoidally and compressively forced simulations 
respectively. The error bar represents the 1-a dispersion. The colored lines represent the 
calibration obtained from each fBm, with density PDF of standard deviation indicated in 
the legend. 
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Fig. 5. — Difference between the PCA slope apcA and the PCA slope afcA predicted by the 
calibration derived from fBms with 0"^ < 2 as a function of log-density dispersion 0"^. The 
trends were obtained from the fBms with lognormal density PDFs and from spectral maps 
generated from rescaled HD density fields combined with fBm velocity fields (see Section 
I6.2p . The shaded area represents the 1, 2, 3o" uncertainty in a^pcA from darkest to lightest. 
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Fig. 6. — One- dimensional ACFs of the first nine principal components (1** to 9*^*) of one 
snapshot of the hydrodynamic simulation with solenoidal forcing (crosses). The solid line 
indicates the best exponential fit to the ACF. The e-fold length (in pixels) is indicated in each 
panel. An exponential ACF describes the ACF of the principal components of a turbulent 
field very accurately for orders > 1. 
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Fig. 7. — PCA pseudo structure function derived from a snapshot PPV cube of hydrody- 
namic simulations (black crosses), from the same simulated cube convolved with a Gaussian 
beam of FWHM 48" (blue triangles), and from the beam-convolved simulated cube after 
scale correctiori (red st ars). For comparison, the scales corrected with the prescription from 
Brunt fc Heyed (J2002al ) are shown as green diamonds. The spatial and velocity scales were 
rescaled to physical units assuming a 22" grid and a temperature of 10 K (sound speed of 
0.2 km s-^. 
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Fig. 8. — Nine first principal components for molecular cloud GRSMC G053. 59+00. 04, 
randomly selected from our sample of 367 molecular clouds from the Galactic Ring Survey 
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Fig. 9. — PGA pseudo structure function for molecular cloud GRSMG G053. 59+00. 04. The 
order of the principal component for each pair of spatial and spectral scales is indicated next 
to each data point. The vertical dashed line shows the resolution limit. Scales detected in 
the 5*^^ and 6**^ are smaller than the resolution limit after scale correction, but above it before 
the correction and thus need to be included in the fit. The solid line represents a power-law 
fit, the slope of which is indicated in the Figure. 
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Fig. 10. — Histograms of the slope of the PCA pseudo structure function obtained from 
GRS clouds, and the exponent P^ of the turbulent spectrum obtained from the calibration 
derived from fBms with purely lognormal PDFs. The errors in the legend correspond to the 
standard deviation of the distributions. The black histogram was derived using the FWHM 
of the beam as the resolution limit (fiducial case). The purple, red, and blue histograms 
show the histogram of apcA derived with resolution limits defined as the la, 3a and 2x 
FWHM widths of the beam respectively. The corresponding mean PCA slopes and /3y are 
also indicated for each case. 
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Fig. 11. — Composite PCA pseudo structure function (composed of all the spatial and 
velocity scales detected in all 367 GRS molecular clouds) shown as a density of points. The 
dashed line indicates a power law of slope 0.62, the average slope of the PCA pseudo structure 
function in the GRS sample, while the solid line shows a bisector fit with slope apcA = 0.6. 
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Fig. 12. — Histogram of the log-density dispersion of GRS m olecular clouds estimated fr om 
Equation [T6| based on their ID velocity dispersion derived in iRoman-Duval et al.l (120101 ). 
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Fig. 13. — Composite PCA pseudo structure function for GRS molecular clouds identified in 
the GRS, with the spatial scales left as angular scales uncorrected for beam smearing. The 
vertical lines indicate several characteristic values of the gaussian beam. 



